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Deep inelastic scattering cross section measurements previously published by the HI and ZEUS collaborations 
are combined. The procedure takes into account the systematic error correlations in a coherent way, leading to a 
significantly reduced overall cross section uncertainty. The combined HI and ZEUS measurements of the inclusive 
neutral and charged current cross sections are used to perform a common NLO QCD fit. The consistent treatment 
of systematic uncertainties in the joint data set ensures that the resulting set of parton density functions (PDFs) 
have a much reduced experimental uncertainty compared to the previous PDF extractions performed separately by 
the HI and ZEUS collaborations. 
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Deep inelastic scattering (DIS) of leptons off nucleons has been essential for our understanding of the structure of 
the nucleon. Colliding electrons and positrons with protons, HERA at DESY has played a unique role in revealing 
details of the proton substructure of quarks and gluons. In June 2007 the operation of HERA was terminated 
after a period of 15 years of data taking. During this time the HI and ZEUS collaborations at HERA successfully 
operated their general purpose detectors which were adapted in 2001 to the luminosity upgrade of the collider 
^ : interaction regions. Both collaborations have published accurate measurements of e^p DIS cross sections based on 
^ , about 120 pb^^ of ep data collected until 2000 [l], 3, H, 0, 0, S 0, IE 0, E3|- These measurements have been crucial 
' in substantially constraining the parton density functions (PDFs) of the proton in global next-to-leading-order 



X 



|-^- (NLO) QCD fits [ill, bj, bj, which use fixed target data as well as HERA data. The HI and ZEUS collabo- 
rations have both demonstrated their capability to extract complete sets of PDFs from their data alone 0, [3] 
[ achieving a precision competitive to that obtained in the global QCD fits. Aiming for the most accurate deter- 
\^ mination of PDFs, HI and ZEUS are now working towards combining their datasets and joining their QCD analyses. 

^ 2. Cross Sections, Structure Functions and Parton Distributions 

The kinematics of lepton hadron scattering is described in terms of the variables Q^, the four-momentum transfer 
' squared of the exchanged vector boson, Bj0rken x, the fraction of the momentum of the incoming nucleon carried by 
I the struck quark (in the quark-parton- model, QPM), and inelasticity y, which is a measure of the energy transfered 
^ ] between the lepton and the nucleon. The differential cross section for the neutral current (NC) process, e^p e^X , 
is given in terms of the structure functions by 

± 

NC _ ^:i::^^Y+F2 T Y^xF^ - i/Fl), (1) 



dxdQ^ ~ xQ^ 



where a is the fine structure constant, 1± = 1 ± (1 — j/)^, and F2, xFs, and Fl denote the generalized NC proton 
structure functions which include Z boson exchange contributions. The structure functions F2, xF^ are directly 
related to quark distributions, and their dependence is predicted by perturbative QCD. The dominant contribution 
to the cross section is F2, which in turn is dominated by the electromagnetic structure function i^2- In the QPM, F2 
- below the bottom mass threshold - can be written as 

4 1 - 

^2 = g [x{u + c + u + c)] + - [x{d + s + d + s)]. (2) 

For low X, X < 10~^, F2 is dominated by the sea quarks, while its evolution is controlled by the gluon contribution. 
Thus HERA data provide crucial information on the low- a; sea-quark and gluon distributions. At high Q^, the 
structure function xF^ becomes increasingly important, and provides information on the valence quark distributions, 
— u — u and dy — d — d. 



1 



34* International Conference on High Energy Physics, Philadelphia, 2008 



The charged current (CC) interactions, e'^{e~)p T'{v)X^ provide additional flavor-type specific information to 
separate quark and anti-quark distributions at high- a; , since their (leading-order) cross sections are given by 

2 

x[{u + c) + {\~yf{d + s)\, (3) 

2 

x[{u + c) + {l-yf(d + ^)\, (4) 

with Gf being the Fermi coupling constant and My/ the mass of the exchanged W boson. 

The e^p NC and CC cross sections can thus be written completely in terms of up-type, xU = x{u + c)^ down-type, 
xD = x{d + s), and of anti-quark type, xU = x{u + c) and xD — x{d + s)), distributions. These (anti-)quark type 
distributions can be determined from HERA data alone. This has the advantage that there is no need for heavy 
target corrections, which must be applied to the i^Fe and fiD fixed target data used in the global fits. Furthermore, 
there is no need to assume isospin symmetry, i.e. that d in the proton is the same as u in the neutron, since 
information on the d distribution can be obtained directly from CC e+p data. 
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3. Combination of H1 and ZEUS Cross Sections 

The HI and ZEUS collaborations have both used their data to perform PDF fits [i, 14|. Both these data sets have 
small statistical errors, so that the contributions of systematic uncertainties become increasingly important and their 
proper treatment is essential. Basically there are two approaches, the Offset method used in the ZEUS analysis and 
the Hessian method applied by the HI analysis. The resulting ZEUS and HI PDFs are compatible, although the 
gluon PDFs do have somewhat different shapes, as is shown in Fig. [T^. 

It is possible to improve on this situation by averaging the HI and ZEUS datasets in a model-independent way 
prior to performing a QCD analysis on them. Since both experiments measure the same cross section at a given 
kinematic point one can minimize the following a-nd obtain a common true value at any measured kinematic point 
together with the correlated systematic shifts 
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Here is the measured central value, and 6i the statistical and uncorrelated systematic uncertainty of the quantity 
M. The M*'*'"'"' are the values following from the minimization; Aaj are parameters for the j*'* source of correlated 
systematic uncertainty and OM'^/daj denotes the sensitivity of point i to source j. For the cross section measurements 
the index i labels a particular measurement at a given (cc, Q^). Equation [5] represents the correlated probability 
distribution function for the quantity M*'*''"^ and for the systematic uncertainty Aaj. 

The defined in Eq. [5] is suitable for measurements in which the uncertainties are absolute or additive, i.e. do 
not depend on the central value of the measurement. For the cross section measurement, however, the correlated and 
uncorrelated systematic errors are proportional to the central values. This proportionality can be approximated by 
a linear dependence. In this case the combination of the data sets using Eq. [5]has a bias towards lower cross section 
values since the measurements with small central values have smaller absolute uncertainties. An improved can 
be defined by replacing Si m' — ^^id^^Ao! — > ^ — which translates the relative or multiplicative 
uncertainties for each measurement to the absolute uncertainty. 

The correlated systematic uncertainties are floated coherently such that each experiment calibrates the other one. 
This allows a significant reduction of the correlated systematic uncertainty for much of the kinematic plane, as is 
shown in Fig. [T)d for three representative a;-bins. 

A study of the global /NDF of the average and of the pull distribution provides a model-independent consistency 
check between the experiments. The HI and ZEUS data are found to be consistent, thus allowing to calculate the 
experimental uncertainties of the PDFs using the tolerance, Ax^ = 1. This represents a further advantage 
compared to those global fits where increased tolerances of Ax^ ^ 1 had to be introduced to ensure that all input 
data sets are consistent with the result of the global fit at a 90% confidence level. 
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HERA I e*p Neutral Cun-ent Scattering - HI and ZEUS 
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Figure 1: Comparison of the PDFs and their uncertainties from ZEUS and HI individual QCD analyses at a fixed value of 
— 10 GeV^ (a) and (b) neutral current e^p cross sections measurements for three selected x bins as function of Q^. The HI 
(open points) and ZEUS data (open squares) are compared to the HI and ZEUS combined data (full points). Measurements 
from the individual experiments have been shifted for clarity. The error bars show the total uncertainty. The curves are NLO 
QCD fits as performed by HI and ZEUS to their own data. 



4. QCD Analysis 

QCD predictions for the structure functions are obtained by solving the DGLAP evolution equations at NLO in the 
MS scheme with the renormalization and factorization scales chosen to be . The DGLAP equations yield the 
PDFs at all values of provided they are input as functions of x at some input scale Qq- For our central fit we 
parameterize the gluon, xg^ the valence quarks, xu.^ and xdv^ as well as the sea anti-quark type xU and xD using 
the generic form 

xfi{x) ^ AiX^'{\- x)'^' [l + DiX + EiX^ + FiX'^) with i = g, u„, d„,I7,U (6) 

and the number of parameters is chosen by saturation of the x^, such that parameters Di, Ei, Fi are varied, i.e. 
they are allowed to take values different from during the fit, but are incorporated only if there is a significant 
improvement of the x^- This leads us to set the parameters Z^i, Ei, Fi = 0, for all partons except xu^ for which only 
Fu^ = 0. The normalization parameters A^^, and Ad^ are constrained to impose the valence quark number sum-rules 
and Ag is constrained to impose the momentum sum-rule. The Bi parameters which determine the \ow-x behavior of 
the PDFs are constrained such that there is a single B parameter for the valance quark distributions -B„^, — Bd^_ = By 
and another B parameter for the sea anti-quark type distributions, B-jj = Bjj = B-q. A flavor decomposition of the 
(anti-) quark type into individual quark flavors can be achieved by assuming that the strange and charm quark at 
the input scale can be expressed as x independent fractions, fs and fc, of I? and U respectively. Imposing that 
d — u— >Oasx^O, a constraint implicit to global fits, further constrains Ajj = Ajy ■ (1 — /s)/(l — fc)- The value 
of fs = 0.33, has been chosen to be consistent with determinations of the strange fraction using neutrino-induced 
di-muon production. The charm fraction, fc = 0.15, has been set to be consistent with dynamic generation of the 
charm from the starting point of = rri^ in a zero-mass- variable-flavor- number scheme. In total 11 parameters of 
the PDFs at an input scale Qo = 4 GeV^ are obtained from the fit. 
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Figure 2: PDFs at = 10 GeV^ from the HI and ZEUS combined data, HERAPDFO.l, compared to recent global PDF 
analyses, CTEQ6.5M (a) and MSTW08 [li (b). 

5. Results 

The set of PDFs obtained from a NLO QCD analysis performed on the combined data set of e^p NC and CC cross 
sections, discussed above, is shown in Fig. [2] (labeled HERAPDFO.l). The much reduced experimental uncertainties 
of the combined cross sections propagate to the error bands of the PDFs, which also include six sources of model 
uncertainties due to variation of mc, rrib, /s, /c, Qqj '^^'^ Qmim the minimum of the data included in the fit. 
The central fit, which in addition to all systematic sources of the combined data set, takes into account 4 sources 
of uncertainties from the combination procedure achieves an excellent x^/NDF of 477/562. Figure [5] also compares 
the HERAPDFO.l to recent global fits. Note that the HERAPDFO.l results employs new, stih preliminary HERA 
data of improved accuracy as discussed above. A further difference to the global fits is that these involve data from 
a much larger variety of experiments and physics processes and thus are forced to use an error definition unlike 
the Ax^ = 1 criterion applicable to the HERA data alone. 

The preliminary PDF set presented here is obtained from the combined, almost complete set of cross section 
measurements of the pre-upgrade data taking period. The current PDF set is available at the LHAPDF library 
and will be published soon, when the HERA I neutral current cross section data are released. The data and QCD 
analyses thus provide a valuable input for the physics at the Tevatron and the LHC. 
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